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ABSTRACT 

Continuum intensity observations obtained with the Michelson Doppler Imager 
(MDI) on-board the SoHO mission provide long time series of filtergrams that are 
ideal for studying the evolution of large-scale phenomena in the solar atmosphere and 
their dependence on solar activity. These filtergrams, however, are not taken in a pure 
continuum spectral band, but are constructed from a proxy, namely a combination of 
filtergrams sampling the Ni I 6768 A line. We studied the sensitivity of this continuum 
proxy to the shape of the nickel line and to the degradation in the instrumental trans- 
mission profiles. We compared continuum intensity measurements in the nearby of 
nickel line with MDI proxy values in three sets of high resolution spectro-polarimetric 
data obtained with the Interferometric Bidimensional Spectrometer (IBIS), and in syn- 
thetic data, obtained from multi-dimensional simulations of magneto-convection and 
one-dimensional atmosphere models. We found that MDI continuum measurements 
require brightness corrections which depend on magnetic field strength, temperature 
and, to a smaller extent, plasma velocity. The correction ranges from 2% to 25% in 
sunspots, and is, on average, less than 2% for other features. The brightness correction 
also varies with position on the disk, with larger variations obtained for sunspots, and 
smaller variations obtained for quiet sun, faculae and micropores. Correction factors 
derived from observations agree with those deduced from the numerical simulations 
when observational effects are taken into account. Finally, we found that the inves- 
tigated potential uncertainties in the transmission characteristics of MDI filters only 
slightly affect the brightness correction to proxy measurements. 

Subject headings: Sun: photosphere - Sun: surface magnetism - Techniques: image 
processing 
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1. Introduction 



The continuum intensity measurements obtained with MDI have provided extensive time se- 
ries of data unaffected by seeing for over 13 years. These data are ideal for studying large-scale 
phenomena in the solar atmosphere and their dependence on solar magnetic activity. They have 
been utiliz ed in many investigat i ons, c oncerning the analysis of large-scale patterns in plasma mo- 
tions (e.g.. iMeunier et al.ll2008L I2007L and references there in), the measurement of the radiative 



prope rties of magnetic elements over the activity cycle ( e.g.lOrtiz et a 



2007 ). and the modeling of irradiance variations (e.g. iKrivova et al 
20091) . 



2002Ll2006l:lMathewetal 



20031 : IWenzler et all 12006 . 



Nevertheless, MDI continuum data are a by-product of the instrument, which was designed 
mainly for Doppler measurements. Specifically, intensities in the continuum are derived by the 
combination of five narrow-band filtergrams, obtained with filters sampling a passband of 94 mA 
FWHM, equally spaced by 75 mA around the Ni I 6768 A mid-photospheric line. The filtergrams 
are labeled Fq. . . F4, where Fq is divided over two bands taken near the continuum, on either side of 
the line, Fi and F4 are centered on the wings of the nickel line, and F2 and F3 are centered around 
its core. Even at Fq the instrument does not sample true continuum. Instead, a proxy continuum- 
intensity filtergram i s constructed by combining the five nominal filtergrams in the following way 
Jscherrer etal.lll995l) : 

Ic = 2Fo -I- /depth/ 2 + /ave , ( 1 ) 

where /ave is the average of the five filtergrams and /depth is the line depth, which is given by: 



/depth = J 2m -F^Y + {F2-F4y 



(2) 



The components of sum Ic theoretically have cancelling systematic errors as a function of solar 
velocity, so that the continuum image is claimed to be free of Doppler induced cross-talk at the 
0.2% level dScherrer et al.ll 19951) . 



The accuracy of the MDI continuum estimate (called the MDI method hereafter) is, however, 
inherently limited by the shape of the spectral line and its sensitivity to both strength and inclination 
of the magnetic field, local thermal stratification, and line-of-sight motions of the observed region. 
In addition, MDI measurements may also suffer from subst antial errors bec ause of uncertainties in 
the actual transmission characteristics of the filters utilized (IWachtenl2008b . 



Several previous studie s have focused on the accuracy of magnetic flux estima t es obtained 



from MDI observations (e.g. Berger & LitesI 2003 : iTran etaL 2005; Ulrich et al. 2009; Wang et al 



20091: bemidov & Balthasai]l2009l:lzhendong et allboiol). as \yell as the accuracy of dynamic 



; mea- 



surements ( Wachter et al.ll2006[Rajaguru et al.ll2o'o7 : lwachte3l2008 ). Results indicate a substantial 
underestimate of magnetic flux and spurious contributions affecting Doppler measurements. This 
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lead to recalibration of magnetograms and dopplergrams series, with the most recent results sug- 
gesting that the calibration of these data could be improved yet furth er. By contrast , the ac curacy 



of continuum intensity measurements has only been investigated by iMathew et alj (120071) . using 
spectral calculations based on one-dimensional atmospheric models. Their results indicate that 
MDI continuum measurements underestimate intensity in sunspots, with an error depending on 
both the magnetic field strength and temperature stratification. 

The objective of this study is to further investigate the uncertainties affecting MDI continuum 
intensity measurements, by considering both high spatial resolution spectro-polarimetric observa- 
tions and numerical simulations. To this aim we defined a brightness correction factor and inves- 
tigated its variation with the physical properties of the analyzed solar features and observational 
conditions, such as spatial and spectral scattered light, finite spatial and spectral resolution and 
position on the solar disk. We also investigated possible effects resulting from degradation of MDI 
filter transmission profiles. 

The paper is organized as follows. We describe the data analyzed and provide some details 
of the reduction applied (Sect. 2). Then, we present the spectral synthesis method and the results 
obtained from our measurements and synthesis computations (Sect. 3). We also investigate sensi- 
tivity of results to filter profiles uncertainties (Sect. 4). We finally discuss the results obtained and 
present our conclusions (Sect. 5). 



2. Spectral data 



To investigate the sensitivity of the MDI continuum measurement proxy we employed pro- 
files of the Ni I 6768 A line obtained from two different sources, namely high-resolution spectro- 
polarimetric observations acquired with th e Interferometric Bidimensional Spectrometer (IBIS 
Cavallinil l2006l : iReardon & Cavallinil l2008b at the NSO/Dunn Solar Telescope, as well as with 
theoretical line synthesis calculations, as described below. 



2.1. Observed spectra 

The observed spectra used in this work were extracted from three IBIS data sets. These will be 
referred to as the Pore, Sunspot, and Faint Plage sets hereafter, alluding to the main solar features 
in their field-of-view (FOV). The Pore set is a spectral scan of the Ni I 6768 A line acquired near 
the disk centre with diffraction limited spatial sampling at 0''083/pixel. The Sunspot and the Faint 
Plage sets are Full-Stokes spectropolarimetric scan of the Ni I 6768 line, acquired at jU = 0.88 
and ji = 0.92, respectively, with a spatial sampling of Of'17/pixel. The IBIS spectral profile at 
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Fig. 1. — Top: subfields of 40 arcsec x 50 arcsec of frames from Sunspot (left), Faint Plage 
(middle) and Pore (right) data sets. Bottom: masks depicting features identified according to 
criteria in Table [T] The colour code is the same as in Fig. |2l 



6768 A has a FWHM of 25 mA (|Reardon & Cavallinill2008h and the spectral sampli ng was about 



30 mA for the three sets. The datasets were r educed following standard procedures (jCauzzi et al 



20081 : IViticchie et al.l 120091: lJudge et al.ll2010h which also correct for the instrumental blue-shift, 
and for instrument and telescope polarization in the case of spectropolarimetric data. The pipeline 
includes a MFBD image reconstruction technique (|Van Noort. Rouppe van der Voort & Lofdahl 
2006b to reduce seeing degradation and possible seeing induced crosstalk. 



Note that all observed line profiles suffer from finite spectral resolution and spectral scattered 
light contamination. This last effect was estimated by com parison of observed average quiet Sun 
profile with Atlas rn e asurem ents by iKurucz et al.l (|1984|) . following the procedure described in 
Cabrera Solana et al.l (|2007L their eq. 1). We found a consistent light contamination level of 12% 
on all the three datasets analyzed. 

To estimate the magnetic flux density along the Line-of-Sight (Bios) in ima ge pixels of the 
Sunsp ot and Faint Plage sets, we applied the Center Of Gravity method (COG) by iRees & Semel 
(|1979l) . This method pr ovides a reasonable estimate of the B^os at the formation height of the line 
core as demonstrated by lUitenbroekl (l2003b for several lines of Fe I. We applied the method to syn- 
thetic sunspot profiles from the 3D MHD simulations (described in Sect. 2.2) taking into account 
both spatial and spectral scattered light and found that the COG in the Ni I line can underestimate 
magnetic flux by up to 28% for B^qs >1.5 kG. The B^qs value deduced for each image pixel 



-5- 



was compensated for LOS projection by applying the relation (|Landi & Landolfill2004l p. 660) 

B = Blos/cosy (3) 
where the inclination angle y derives from 



Y = arccos , (4) 

X 

x = V / \/U^ ^ Q^, and x is computed at the left wing of the line. 

As an aid in classification of the observed pixels (see Table [B, we computed the total polar- 
ization signal T in each image pixel, which is defined as 



^ = J{Q^ + U^ + V^y/^/IcdX 



(5) 



where Ic is the intensity in the continuum near the Ni I line. 



Finally, line-of-sight velocities were estimated by Doppler shifts of line cores with respect to 
the average quiet Sun line center position. 

Characteristic maps of the continuum intensity Ic of each set are shown in Fig. [H The FOVs 
include a large variety of solar features, characterized by different physical properties and therefore 
different shapes of line profiles. An example of such variety is reported in Fig. [2l which shows 
average line profiles of pixels selected according to the criteria described in Table [T] 



2.2. Synthetic spectra 



To better understand the results obtained from observations and to distinguish between obser- 
vational and methodological effects, we also computed synthetic spectra of the Ni I line by utilizing 
various numerical atmospheres: two two-dimensional cross sections from a three-dim ensional sim- 
ulatio n of magneto-convection of an active region with average vertical field of 250 G (IStein & Nordlund 
1998 ). a vertical cross section of an umbral magneto-convection simulation ( Vitas et aD 2010 ). and 
various one-dimensional hydrostatic models. The one-dimensio nal static models we considered 
ra nged from the semi-e r npiric al umbral models E, L, and M from lMaltby et al.l (|l986l) . the model 
by lSocas-Navarro et al.l (120041 Socas, hereafter) to the most recent radiative equilibrium models of 
KuruczQ. In all the Maltby and Kurucz models, we assumed vertical height-independent magnetic 
fields. For the Maltby models, we assumed magnetic field strengths ranging from 2.0 kG to 4.0 



http://wwwuser.oat.ts.astro.it/castelli/grids.html 
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Table 1 : Criteria utilized for the identification of the solar features depicted by masks in Fig. [Hand 
whose line profiles are shown in Fig. [2l 



Solar feature 


Criterion 


Data set 


quiet Sun 


L ^ j.yjinA 


Spot 






Faint Plage 


facula 


Blos > 0.2 kG 


Faint Plage 




and Ccore > 10% 




umbra 


Blos> 1.5 kG 


Spot 


penumbra blue sh. 


0.8kG<5LO5< 1-1 kG 


Spot 




and vlos < 




penumbra red sh. 


0.8kG <5lo5 < 1-1 kG 


Spot 




or VLOS > 




pore 


/ < ave — 5(7 


Pore 



Note. — Bios indicates the magnetic flux strength along the line of sight, Ccore is the intensity contrast in the Ni I 
line core that is defined by the ratio I/Iqs, where Iqs is the quiet Sun intensity, vios indicate the velocity of plasma 
along the LOS, ave and a are the average and standard deviation of all the intensity values measured on the FOV. 




Fig. 2. — Line profiles, computed by Gaussian fits to average spectral intensities of various fea- 
tures singled out on the IBIS observations as indicated in Table 1 . Continuous Black: quiet Sun. 
Continuous Magenta: facula. Dashed Green: umbra. Dot-dashed Blue: pore. Triple-dot dashed 
Light blue: blue shifted penumbra. Triple-dot dashed Red: red shifted penumbra. 
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Fig. 3. — Line profiles derived from the numerical simulations. Top panel: results from the active 
region magneto-convection simulation. Solid: quiet Sun. Dashed: upflow. Dot-dashed: bright 
magnetic feature (facula). Dot dot dashed: dark magnetic feature (micropore). Bottom panel: 
results from hydrostatic models and magneto-convection simulations of sunspot. Overlaied profiles 
show MDI spectral transmission profiles normalized to fit the figure 
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kG. The full set of one-dimensional models, which are all representative of umbral regions, and the 
corresponding magnetic field strengths are listed in Table [21 The field strength of the Socas model 
listed in this table, however, is the value at t = 1 at 5000 A, since this model was used with its 
original depth-dependent magnetic field structure. Note that we set the turbulent velocity value in 
Kurucz models to zero, since other values in these models are more representative of non-magnetic 
and convectively unstable atmospheres. 

Th e spectral synthesis was p erformed in non-local thermal equi librium (NLT E) with the RH 
code by lUitenbroekl (120021 l2003h . and employing the model atom of lBrulsl (|1993[) . The emergent 
spectra were calculated at various lines of sight for a 1 .A wide spectral interval centered on the 
Ni I line. Figure [3] shows examples of line profiles derived from the synthesis in th e various atmo 



sphere s, all normalized to the intensity value calculated from the FAL-C model of iFontenla et al. 



(Il993b . which is representative of the quiet Sun. This normalization provides a common intensity 
reference value for the different simulations. Similarly to observations, the line profiles derived 
from the synthesis display a large variety of shapes, because of the different physical conditions 
in the simulated atmospheres. The synthetic profiles differ from those deduced from observations 
of magnetic features mostly in ter ms of the l ine core reversal. This must be mainly ascribed to 
the transverse Zeeman component (jStix 1120041) . which, especially for high values of magnetic field 
strength, can considerably modify the shape of the line core, and the fact that all observations were 
taken slightly off disk center. 



3. Simulation of MDI continuum measurements 

To test the algorithm for the MDI continuum proxy we compare the values it predicts with the 
real continuum values in the observed and simulated spectra of the Ni I 6768 line. We first evaluate 
the intensities in each of the five filter positions by multiplying the spectra with the filter passband 
and integrating over wavelength, then add the integrated intensities according to Eq. [TJ The MDI 
transmission profiles were centered on the core of the average line profile of the Quiet Sun. This 
was defined as in Table 1 in the case of Sunspot and Faint Plage dataset, while for the Pore dataset 
we considered an area in the FOV of 80 square arcsec far from the pores. As for simulations, we 
centered the transmission profiles to the core of the line profile derived from the FAL-C model. 

Since we generally do not measure absolute intensities in observations, but rather intensities 
relative to some average intensity, we analyze validity of the MDI continuum estimate by measur- 
ing contrast, where contrast is defined as the ratio between the intensity of an analyzed feature and 
the intensity of the quiet Sun. For observed spectra, the contrast was estimated after compensation 
for the center-to-limb variation (CLV) of quiet Sun intensity across the FOV on each frame. The 
CLV was derived from the second order polynomial surface best fitting the intensity pattern of quiet 
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Table 2: One-dimensional atmosphere models considered. Kurucz models are indicated with 
Kur_XXX where XXX is the Effective temperature in 10-^ K units. 



Model B(kG) 



Kur_5.50 


1 


.10 


Kur_5.25 


1 


.70 


Kur_5.00 


2, 


.25 


Kur_4.75 


2, 


.25 


Kur_4.50 


2, 


.50 


Kur_4.25 


3, 


.00 


Kur_4.00 


3, 


.50 


Kur_3.75 


4, 


.00 


Maltby L 


2, 


.00 


Maltby L 


2, 


.50 


Maltby L 


3, 


.00 


Maltby M 


2, 


.50 


Maltby M 


3, 


.00 


Maltby M 


3, 


.50 


Maltby E 


3, 


.00 


Maltby E 


3, 


.50 


Maltby E 


4, 


.00 


Socas 


2, 


.3 



Sun regions in the image, after masking dark magnetic features by applying an intensity criterion 
thresholding, and application of a smoothing function to pixel intensities. Concerning numerical 
simulations, intensity values were normalized by the intensity at the Ni I continuum calculated 
with the FAL-C model for the appropriate disk positions. Furthermore, we define as brightness 
correction factor to MDI continua measurements the ratio C/Cmdu where C is the contrast value 
obtained in the continuum near the Ni I line (hereafter called the continuum contrast) and Cmdi is 
the contrast derived with the MDI method (hereafter called the MDI contrast). In the following we 
investigate the dependence of this correction factor on various properties of the analysed regions, 
e.g. the contrast, magnetic field strength, and plasma velocity, as well as observational conditions. 



3.1. Results from observed spectra 

Figure m shows maps of the intensity contrast {left panels) in the Sunspot {top) and Pore {bot- 
tom) sets and the corresponding maps of the brightness correction to MDI measurements {right 
panels). We found that the correction factor is, on average, lower than 2% of the MDI contrast, 
but increases up to 8-10% in regions with high magnetic field strength. The figure also shows that 
the MDI contrast is, on average, overestimated in upflow regions (granules), while it is underes- 
timated in regions with downflows (intergranular lanes), and in particular in regions with strong 
magnetic fields. The map of the brightness correction in the Pore {lower right panel) looks very 
similar to a map of the inverse granulation as it would appear in core intensity of a photospheric 
spectral line, implying that the correction factor depends strongly on the photospheric temperature 
gradient, which is steep over granules and shallow over inter- granular lanes. The reason that the 
MDI continuum estimate depends on this temperature gradient lies in the corresponding widths of 
the spectral line profile, which are narrow and deep in granules, and wide and shallow over inter- 
granular lanes. Profiles that are wider than that of the average quiet Sun cause an underestimate of 

0.33 0-7-1 1,15 0.98 ' i.O.? lO.T^, ' 1.10 0.99 1.0' 1.02 




Fig. 4. — Left panels: continuum contrast maps derived from Sunspot (top) and Pore (bottom) 
datasets. Right panels: corresponding maps of the brightness correction. 
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Fig. 5. — Dependence of the brightness correction factor on the continuum contrast derived with 
MDI method. Red asterisks, blue triangles, and green diamonds show values derived from averages 
over contrast bins, each 0.02 wide, on the Sunspot, Pore, and Faint Plage sets, respectively. The 
error bars represent the standard deviations of the values in each bin. Shaded area delimits the 
values within 3% of the maximum of the distribution of values in each bin for the Sunspot data set. 
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Fig. 6. — Dependence of the brightness correction on the magnetic flux of the analyzed region. 
Legend as in Fig 151 
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the line depth (Eq. [2l), and consequentially an underestimate of the MDI continuum (Eq. [T]). The 
effect is also visually enhanced by the fact that magnetic field concentrations preferentially reside 
in downflow regions. This is further explained in Section 3.2 below. 

Figure |5] shows the brightness correction to the MDI contrasts derived from all the image 
pixels in all three data sets. Symbols and error bars of a given colour show the average and standard 
deviation of the correction factor derived from a given data set by considering the MDI contrast 
values for the set collected into bins with 0.02 width. We found that the brightness correction is 
close to unity for bright features, because of the opposite contribution from Doppler shifts and 
magnetic fields to the line shape, though we note some dispersion of results. On the other hand, 
the correction factor increases with the decrease of the MDI contrast derived from the images, 
reaching 1.08 for MDI contrasts < 0.3. 

Next we investigated the effects of the magnetic field strength on MDI contrast estimates in 
more detail. Figure [6] shows the results derived from a sample of Sunspot and Faint Plage images 
(7 and 5 IBIS line scans, respectively). We found that the brightness correction increases with the 
increase of the magnetic field strength, although a saturation value of »i 1 .04 is reached at magnetic 
flux larger than 1 .5 kG. In spite of the small standard deviation values, we also notice a large scatter 
of results, especially for the regions with \^{B) | < 1 .5 kG. In the range 0.7 kG < \^{B) \< 1 .2 kG 
this is due to combined effect of LOS velocity and magnetic field strength, since this range also 
comprises pixels belonging to the sunspot penumbra. At lower values of magnetic flux the large 
scatter of results likely derives from uncertainties in the COG estimations of magnetic field, which 
are due to limited spatial resolution and residual cross-talk in Stokes signals. 

Finally, we investigated the effects of LOS velocity on brightness correction, restricting the 
analysis to finite ranges of magnetic field strength values. Results obtained for quiet Sun regions 
selected on the three datasets are shown in Fig. |71 In agreement with visual inspection of images in 
Fig. IHthe brightness correction is, on average, larger than unity for downflow regions and smaller 
than unity for upflow regions, but with large dispersion. Due to the lower image quality of Pore 
and Faint Plage datasets, this effect is more evident on the Sunspot data. 



3.2. Results from synthetic spectra 

Our analysis of observations performed with the observed Ni I spectra indicates that the bright- 
ness correction to MDI measurements depends on the solar feature analysed. This is confirmed by 
the results derived from various numerical simulations. Similarly to Fig. [5] derived from obser- 
vations. Fig. [8] shows the brightness correction deduced from the various atmospheres described 
in Sect. 2. Here the set of points with contrast close to unity stems from the magneto-convection 
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-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 
Vlos (km/s) 



Fig. 7. — Dependence of the brightness correction on the LOS velocity for quiet Sun regions 
identified on the three IBIS dataset. Negative velocities correspond to upflow motions, positive to 
downflow motions. Legend as in Fig. [5] 
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Fig. 8. — Dependence of the brightness correction to the MDI contrast computed for the various 
atmospheres considered in this study. Crosses: results from the magneto-convection calculations. 
Asterisks: Kurucz models. Triangles: Maltby and Socas models. Black: results obtained at disk 
center. Red: results obtained at disk position /i=0.88 after data degradation to account for the finite 
spectral resolution, spurious light and scattered light degradations affecting the IBIS data. 
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simulation of active region, the set of points with contrast close to 0.1 stems from the magneto- 
convection simulation of sunspot umbra and the intermediate contrast points stem from hydro- 
static models. We found that the brightness correction is close to unity, with some dispersion of 
results due to magnetic field and velocity effects, for MDI contrast values ranging from 0.8 to 
1.2, as also derived from the three IBIS data sets. However, the correction factor increases up to 
1.15-1.25 for MDI contrast <0.5, i.e., by a larger fraction than derived from observations. We also 
found a large dispersion of results from the magneto-convection simulations of umbral regions, 
because of the wide range of field strengths in the analysed atmosphere models. The dispersion of 
correction values from simulations is also larger than those from observations. 

To assess the importance of observational effects on the results derived from the three data 
sets and to place them in context with the results deduced from model calculations, we computed 
synthetic spectra from the various model atmospheres at the heliocentric angles of our observations 
(e.g. /i = 0.88 for the Sunspot set) and degraded the results to account for the finite spectral 
resolution, spurious spectral light contamination, and spatially scattered light degradation affecting 
the observations. 

This latter effect was mimicked by applying the formula /5 = (1 — alq, where a in- 

dicates the level of spurious spatial light contamination, and Iq is the quiet Sun intensity value. 
Estimation and compensation on high spatial resolution data for scattered light results uncom- 
plete, even eniployin g the most sophisticated restoration algorithm available (e.g. MOMFBD 



Scharmer et al.l 120101) . In the lack of measure ments, we as s umed the value a = 0.2, which is 



within the range of values (0.17-0.25) found by lBriand et all (120061) for IPM/THEMIS data. The 



results obtained indicate that the observational conditions and degradations above listed all af- 
fect the brightness correction, specifically by decreasing its value. The difference between the 
three IBIS datasets must therefore be ascribed to the different image qualities. We also note that, 
among the observational effects considered, the brightness correction is affected the most by spu- 
rious spatial l ight contamination, which also contributes to a smaller dispersion of the results from 



observations (IMathew et al.ll2007l : ICriscuoli & Ermollill2008b 



Figure ID shows the dependence of the brightness correction on the magnetic field strength in 
the synthetic atmospheres. The results deduced with simulations are qualitatively in agreement 
with those derived from observations, though the saturation value is reached only at B = 3 kG and 
the dispersion of results is larger than obtained from observations as a result of image degradation 
by limited spatial resolution and spatially spurious light. The plot clearly shows a double distri- 
bution: one for which the brightness correction value is smaller than 1 .05 and one for which it is 
higher. We found that data from the active region magneto-convection simulation belong to the 
first group, whereas the second group is populated by sunspot models (with the exception of the 
two hottest Kurucz models) and the umbral simulation. 
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We found that the MDI continuum estimate method depends in a complicated way on both 
temperature and magnetic field intensity. The larger the magnetic field strength, the wider the 
line profile. If the line profile becomes so wide and/or shifted that it reaches into (one or two 
of) the outermost transmission Fq passbands, the continuum (Eq. [B is und erestimated, resulting 



in a brightness correction larger than one. This effect is also discussed in iMathew et al.l (120071) 
and illustrated in their Fig. 7. These authors also notice that for the largest magnetic field values 
considered in their study, a saturation of the error occurs because of the filling in of the line core, 
as is also clear in the line profiles in Fi g. [31 By contrast, an increase in the effective temperature 



narrows the line and increases its depth (IMathew et al.ll2007L their figure 8). As a consequence, for 



a given value of magnetic field strength the brightness correction is smaller for faculae rather than 
for sunspots. On the other hand, at the saturation limit, where the field mostly affects the width of 
the line, the error in the MDI continuum estimate increases with the effective temperature as the 
line narrows. We also note that the brightness correction factor for the synthesized spectra (Fig.|9l) 
seems less affected by saturation than the one for the observed spectra (Fig. (6]). 

Next, we investigated the dependence of the brightness correction on the LOS velocity of the 
plasma. We found that, on average, the brightness correction is larger than unity and increases with 
the magnitude of the LOS velocity. This is illustrated in Fig[10l which shows the results derived 
for quiet-Sun portions (defined here by those pixels with 5 < 50 G) of the magneto-convection 
simulations representative of active region. As we can see, d ata from simulation s allow to investi- 



gate a wider range of velocities with respect to observations (jCauzzi et al.ll2006r) . When restricting 
the analyses to the range from -2 km/s to 2 km/s we note the same asymmetry we found for ob- 
servations, namely a brightness correction smaller than unity for upflows and larger than unity for 
downfiows. This asymmetry is the result of the correlation between temperature gradient and up- 
or downfiows in the granulation. Granules correspond to upflows (negative vlos) and have steep 
temperature gradients resulting in narrow and deep profiles. The shallow temperature gradient over 
inter- granular lanes, where most downfiows occur, leads to broader profiles. The spacing of the 
MDI filter positions is chosen such that it matches the width of the average quiet-Sun profile of the 
Nickel line in the sense that Eq. [2] gives the proper line depth. For substantially narrower profiles 
the filter positions Fi and F4 sample close to continuum values increasing the line depth estimate, 
overestimating the continuum. This leads to the correction factors below unity for small upflows. 
For larger line shifts to either sides, but already for smaller velocities in the downflows with their 
wider profiles, the correction factor is larger than one, because Fi and F2 (for blue shifts), or F3 and 
F4 (for red shifts) both end up close to the core, minimising their difference in sampled intensity, 
underestimating the line depth, and underestimating the continuum level. 

Finally, we considered the dependence of the brightness correction on the position of the 
analysed feature on the solar disk. To this aim, we computed the synthetic spectra for 13 different 
viewing angles from various model atmospheres. In particular, we employed all the static sunspot 
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Fig. 9. — Dependence of the brightness correction on the magnetic field strength derived from 
synthetic spectra. Legend as in Fig. [8] 
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Fig. 10. — Dependence of the brightness correction to MDI measurements on the LOS velocity 
for quiet regions from magneto-convection simulation of active region. 
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models as well as the magneto-convection simulations representative of an active region. The 
Nickel line forms low enough in these models so that the vertical domain includes its formation 
height, even at high inclinations of the LOS. Figure [TT] shows the correction factor for the different 
Kurucz models with cooler models and stronger fields (both leading to wider profiles, and larger 
correction factors) towards the top. For the models with B <3 kG (or Teff > 4500 K), the bright- 
ness correction monotonically decreases with the decrease of the heliocentric angle, because the 
width of the line decreases with the line-of- sight component of the magnetic field, which changes 
linearly with /i for the assumed vertical fields. For models with B > 3 kG ( Tgff < 4250 K), 
the CLV of the correction factor increases strongly with decreasing ji from 1 .0 to 0.9, and then 
decreases again monotonically with decreasing /i. The initial rise is the result of the transverse 
magnetic Zeeman effect, which produces intensity decrease of the core of the profile for positions 
close to disk center, thus creating characteristic double-peaked profiles. These quickly disappear 
for /i < 0.9. Figure [TT] also shows the CLV of the correction factor for faculae (squares), desig- 
nated as those areas that have a field strength larger than 0.2 kG and contrast greater than one) in 
the magneto-convection active region simulation. These facular regions show very little variation 
in their required correction factors accross the disk. 



4. Sensitivity to passband uncertainties 

To assess the importance of an inaccurate knowledge of the filter characteristics, we analysed 
the dependence of the brightness correction derived from observed spectra on potential inaccura- 
cies in the precise tuning and width of the MDI passbands. In particular, we computed the variation 
of the brightness correction for displacements of the Fq.. .F4 passband pattern ranging from -0.14 
A to -1-0.14 A from its nominal position. These values correspond to Doppler shifts of ± 6 km/s 
and are 1.5 times la rger than the ones considered in the calibration of the MDI measurements 



dScherrer et al.lll995L their Fig. 12). We also investigated the sensitivity of the required brightness 



correction to variations of the width of the MDI passband respect to the nominal value of 94 mA. 

Figure [111 displays the resulting brightness correction as a function of displacement of the 
MDI pass band pattern, when considering the same average observed line profiles for quiet Sun, 
plage, pore, umbra, and the most blue-shifted penumbra (see Fig. [2]). We found that, for moderate 
displacement values, the brightness correction decreases for displacements in the same direction 
as the line profile Doppler shift, while it increases for shifts in the opposite direction. This is 
particularly clear for the penumbral profile (triple-dot-dashed curve). For large displacements all 
the curves in Fig. [121 show a trend towards unity, because the offset between the line center and 
transmission pattern becomes so large that the passbands sample more a nd more t r ue con tinuum. 



Note also the oscillatory variations of the curves, that, as pointed out by IWachten (|2008|) . have a 
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period proportional to the spacing of the transmission profiles of the filters. 

Next, we computed variation of the brightness correction for different assumed widths of 
the tunable MDI passband by convolving the nominal profiles with Gaussians of widths that are 
multiples of the width of the profiles. This last quantity was estimated fitting the highest lobe of 
Fq with a Gaussian function. Figure [T2l shows the results for the various solar features. We found 
that the brightness correction decreases with the increase of the filter profile widths, and eventually 
drops below unity for the pore and facular profiles. This is the result of the transmission profile 
becoming wider than the line profile (see also the discussion in 3.2). Note also that the brightness 
correction for the penumbra is the least affected by the widening of transmission line profiles, since 
its value is mostly determined by the shift of the spectral line. 



5. Discussion and Conclusions 

We have presented a study of the accuracy of the MDI continuum intensity measurement, with 
the help of high spatial resolution spectro-polarimetric observations acquired with the IBIS/DST 
and outcomes from numerical simulations representative of different solar regions. In order to 
provide the reader with a quantity that clearly indicates the MDI uncertainties and allows their 
compensation on measurements, we defined as brightness correction the ratio between the contrast 
value in the continuum and the one estimated through MDI method and investigated the variation 
of this quantity with the physical properties of analyzed solar regions. 

We found that the MDI brightness correction, on average, increases with the increase of the 
LOS magnetic field strength, the decrease of temperature, and the increase of LOS velocity. To 
summarize, the contrast is overestimated in quiet upflow regions, and is underestimated in mag- 
netic and quiet downflow regions. The amount of correction is less than 1% in quiet regions and 
increases with the magnetic field strengths, reaching values up to 25% of the MDI continuum value 
in sunspots. However, line saturation effects limit the MDI error in regions with highest magnetic 
field strengths. 

Results derived from observations agree with those deduced from simulations, when ac- 
counting for the observational conditions and instrumental degradations affecting the observations. 
Among the effects considered, spatially spurious light is the one that affects the most the bright- 
ness correction values. In the lack of a reliable estimate on our observational data, we therefore 
refrained from a direct comparison of results from MDI and IBIS data. On the other hand, we did 
not find any significant variation of the brightness correction values when degrading the results 
from IBIS data to the spatial resolution of MDI, thus suggesting that brightness correction is not 
significantly dependent on the spatial resolution of the observations. 
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Note that the brightness correction values that we found correspond to errors much higher 



than the 0.2% quoted in the instrument description of IScherrer et alj (119951 p. 163). 



The brightness correction factor derived from synthetic spectra also depends on the position 
on the solar disk in different manners for the different features analyzed. We found that the bright- 
ness correction for sunspot models with strong field (5los ^ 3 kG) varies by up 2% for disk 
positions corresponding to 0.9 < /i < 1, reaching a maximum value for /i = 0.9, and then decreas- 
ing towards the limb. For weaker fields the variation across the disk is much less (Fig. [TT]) . Also 
for quiet Sun, faculae and micropores the variation in brightness correction is small, below 1%, 
across the disk. 

Finally, we investigated effects of possible filter degradations on MDI continuum contrast 
estimates. We estimated variations with the shift of the passband pattern and with the change in 
the width of the transmission profiles of the instrument. Results are summarized in Fig. [12] for 
various features observed on the solar photosphere. We found that for small amounts ([-0.03,0.03] 
A) of transmission profiles displacement the variation of the brightness correction is negligible for 
features not affected by Doppler shifts (i.e. faculae and quiet sun). For line profiles shifted with 
respect to quiet sun line center, a shift of filter center in the same direction as the line Doppler shift 
causes a reduction of the amount of correction, whereas a displacement in the opposite direction 
causes an increase. In addition, we found that broadening of the passband causes a decrease of 
the amount of correction when the width of the line profiles is larger than that of the passband. 
For transmission line profiles wider than the profile of the solar feature under study, the contrast is 
overestimated and the amount of deviation of the correction coefficient from unity increases with 
the increase of the filter degradation. 



Our results are qualitatively in agreement with those pres ented b y Math ew et al.l (l2007h . who 



applied MDI procedure to synthetic line profiles obtained from lKurucz I (i 1 99 ll) atmosphere models. 
Nevertheless, a comparison of results shows that they obtained brightness correction larger than the 
ones we have presented. These discrepancies must be ascribed to the differences in model atmo- 
spheres employed, assumptions for atomic parameters and micro and macro velocity values. Due 
to the arbitrariness of the choice of these quantities inherent to ID models, this analyses should 
be repeated by a direct comparison of MDI data with observations in the Ni I continuum prop- 
erly compensated for scattered light, or with the employment of magneto-convection simulations 
encompassing a larger variety of sunspots. 

Results presented in this study show that photometric measurements obtained with MDI data 
should be revised by taking into account the brightness correction factors presented. Due to the 
dependence of this factor on observational conditions, continuum images should be compensated 
for spatial scattered light before any correction is applied. Since MDI Point Spread Function varies 
with time, this result contributes to question studies based on the temporal variation of photometric 



properties of solar features derived from MDI data analyses. 

Finally, our study proofs forward modeling as a powerful tool also for instrumental calibra- 
tion. In particular, results presented here are of interest for the interpretation of data acquired with 
the Helioseismic and Magnetic Imager onboard the Solar Dynamic Observatory. 
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Fig. 1 1 . — Center to Limb variation of the brightness correction for various magnetic features 
investigated through synthetic spectra. Squares: faculae, from magneto-convection simulation. 
Asterisks: sunspots, from Kurucz models. 
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Fig. 12. — Variation of brightness correction of features identified on IBIS data (see Table [T]) 
with the filter drifts (top panel) and modifications of shapes of transmission profiles (bottom panel) 
considered in this study. Continuous line: plage region. Dashed line: umbra. Dot-dashed line: 
pore. Triple-dot-dashed line: penumbra. 



